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Abstract 

We derive effective equations of motion for a massless charged particle coupled to the 
, dynamical electromagnetic field having regard to the radiation back reaction. It is shown that 

OO , unlike the massive case not all the divergences resulting from the self-action of the particle are 

, Lagrangian, i.e. can be canceled out by adding appropriate counterterms to the original action. 

, Besides, the order of renormalized differential equations governing the effective dynamics turns 

. out to be greater than the order of the corresponding Lorentz-Dirac equation for a massive 

particle. For the case of homogeneous external field the first radiative correction to the Lorentz 
equation is explicitly derived via the reduction of order procedure. 



1 Introduction 



Oh! 

X 

■ One of the interesting features of the classical electrodynamics is a possibility to account for the 
self-action of a point charge in a purely local manner by adding certain higher-derivative terms to 
the usual Lorentz's equations. Physically these terms are responsible for the back reaction of the 
radiation emitted by an accelerating charge. The first systematic derivation of the effective equations 
of motion for d = 4 relativistic massive particle was given by Dirac Since the non-relativistic 
limit of Dirac's equation coincides with the zero-size limit of Lorentz's model of an electron [0 
it is often referred to at the Lorentz-Dirac equation. The case of d = 4 electrodynamics is not 
an exceptional one and the analogous effective equations can be obtained for a massive particle in 
arbitrary dimensions, but the local equations take place only for even d^s. The case of d = Q massive 
particle was first considered in Ref. The general framework to the Lorentz-Dirac equation in 
higher dimensions was developed in our recent paper 0]. 

Although there is a large literature devoted to different aspects of the Lorentz-Dirac equation, 
its solutions and applications (see p, 0, § for a modern review and further references) the case of 
a massless particle has somehow escaped consideration. The aim of this paper is to fill this gap in 
the case of d = 4 massless particle. We emphasize that the effective dynamics of massless particle 



is essentially different from that of the massive one, and it can't be obtained as a massless limit of 
the Lorentz-Dirac equation. 

While massless electrically charged particles (so to speak, massless electrons or charged photons) 
have not been yet observed experimentally no obvious prohibitions are known against their existence 
at least from the viewpoint of classical field theory. Moreover, these particles are predicted by 
supersymmetric gauge theories with an unbroken sector of Abelian gauge symmetry. Perhaps some 
peculiarities of the effective dynamics, revealed in this paper, do indicate a certain inconsistency of 
the classical electrodynamics of massless particles with radiation effects included. 

This paper may also be viewed as a preparative for study the radiation reaction problem in the 
case of extended objects like strings and branes universally coupled to the dynamical antisymmetric 
tensor fields. It is notable that the regularization procedure we use here is closely analogous to the 
formalism of characteristic (de Rham's) currents intensively exploited for the study of anomalies in 
5-branes coupled to 11-dimensional supergravity background |12, 13, |l4i . 

The paper is organized as follows. In Sec. 2 we discuss some peculiarities of the dynamics 
of a relativistic massless particle coupled to external electromagnetic field and describe a class of 
isotropic word lines for which the radiation reaction effects can be properly treated. Taking into 
account the radiation back reaction leads to inevitable infinities resulting from the "pointness" of 
the particle. The renormalization procedure for these classical infinities is briefly explained in Sec. 
3. Here we also present our main result - the effective equations of motion for a massless charged 
particle, and trace the origin of non-Lagrangian divergences. As for the Lorentz-Dirac equation, 
our equations involve higher derivatives and thus the reduction-of-order procedure is required to 
assign them with a proper mechanical interpretation. This is done in Sec. 4. In particular, we 
explicitly derive the leading correction to the usual Lorentz's equation for the case of homogeneous 
external fields. In the concluding section we summarize the results. Appendix contains details of 
the calculations omitted in Sec. 3. 



2 Massless particle in the classical electrodynamics 

Let M^'^ be four-dimensional Minkowski space with coordinates {x^}, /i = 0, 1, 2, 3, and the metric 
Vfiu = diag{—l, +1, +1, +1). Consider a massless point particle moving in M^'^ and coupled to the 
electromagnetic field. The dynamics of the whole system (field) + (particle) is governed by the usual 
action functional] 

S = -\j drex^ + J d'x (^A^j ^ - -^F^^F^^^ . (1) 

Here F^^, = d^A^, — d^A^ is the strength tensor of the electromagnetic field, e 7^ is the Lagrange 
multiplier and 

f{x) = q j 5^{x - x{T))x^'{T)dT (2) 
^We take unified dimensions for space and time setting c = 1. 



is the density of the electric current produced by the point charge q moving along a world line 
x^{t) with the four- velocity = dx^/dr. The action is invariant under arbitrary world-line 
reparametrizations t ^ t' , provided e(r) is transformed as a world-line density (einbein) 

e(r') = ^e(r) . (3) 

Varying the action (|l|) with respect to e one gets the standard isotropy condition (Lagrangian 
constraint) on the four-velocity of the massless particle 

x^ = 0, (4) 

i.e. a tangent vector to each point of particle's trajectory x^{r) lies on the future light cone attached 
at this point. We will refer to such a trajectory as an isotropic one. The variations of the action 
with respect to particle's trajectory x^{t) and electromagnetic potentials A^{x) lead to the coupled 
systems of the Lorentz and Maxwell equations 

ex^ + ex^ = qF^^x", (5) 



dT^, = Anj^ . (6) 

Using the reparametrization invariance of the model and the transformation property one can 
always bring the einbein e into the form e = 1 or, equivalently, one may impose this condition to 
fix a particular parametrization of the world line. In this gauge the form of the Lorentz equations 
is quite similar to that for a massive particle, 

•^/i qF^iiX , (7) 

written in a proper-time parametrization (i;^ = 1). The crucial difference, however, is a presence 
of the additional isotropy condition (^) having no analog in the massive case. In particular this 
suggests impossibility to introduce the proper-time parametrization for world lines of the massless 
particle. 

Another convenient parametrization to use when working with the isotropic world lines is ex- 
tracted by the following gauge fixing condition: 

x^ = l. (8) 

One may regard it as a "massless" counterpart of the proper-time parametrization. Contrary 
to the massive case, however, there are some isotropic trajectories which do not admit such a 
parametrization. Since for any isotropic curve > 0, the trajectories for which = at some r's 
may be thought of as degenerate ones in the sense that the complementary set of nondegenerate 
trajectories fills an open and everywhere dense domain in the space of all isotropic curves, so that 
any degenerate isotropic curve can be made nondegenerate by a small perturbation. 

o 



On the other hand, if for some interval Ti < r < T2 the four-vector of acceleration is know to be 
isotropic, = 0, then we can immediately conclude that this part of the trajectory is given by a 
straight segment. Indeed, since the velocity and acceleration of the particle are always orthogonal 
to each other we get a pair of isotropic and orthogonal vectors x'^ and x^, but any such two vectors 
are known to be proportional to each other in the Minkowski space, 

x^' = X{t)x''. (9) 

Integrating the last equation we see that the particle does move along a straight line until ti < r < 

The positive-curvature condition (x^ > 0) for particle's trajectory implies a nonzero value of 
the external electromagnetic field in the right-hand side of the Lorentz equation @, otherwise the 
particle would move along a straight isotropic line for which x^ = 0. More precisely, squaring both 
sides of the Lorentz equation and accounting that, by definition, e 7^ we arrive at the following 
condition: 

T^,(x)i;^x" > 0, (10) 

where 

T,u = (f/F,. + \v,uF''^F^P^ (11) 

is the stress-energy tensor of the electromagnetic field. Passing to an appropriate Lorentz frame it 
is not hard to see that T^^x^x^ > for any field F^^ = (E, H), and the equality implies either 

V II E II H (12) 

or 

v||ExH and E ■ H = , E^ = H^ (13) 

V = {x*} being the space part of the four-velocity. The latter situation may occur, for example, for 
a massless particle moving on a background of free electromagnetic wave. 

In this paper we restrict our attention to the self-consistent solutions for the Lorentz-Maxwell 
equations (^, |]) obeying the gauge fixing condition (||) and the compatibility condition (|1^). In the 
next section we will see how these restrictions arise naturally upon calculating the radiation-reaction 
force. 

3 Radiation reaction and renormalization 

In this section we follow the approach of work |^ to account for the radiation back reaction to 
the relativistic motion of a point charge. For this end, we first solve the Maxwell equations (|^) 
expressing the electromagnetic field as a functional of particle's trajectory. The effective equations 
of motion are then obtained by substituting this solution to the Lorentz equations (^). Inasmuch 
as the electromagnetic field is singular at the points of particle's trajectory the last step requires a 
renormalization procedure to remove inevitable infinities. The alternative approach to the problem 



(also involving some renormalization) is based on the computation of energy radiated by an accel- 
erating charge and identification of this energy with the work of a hypothetical (radiation reaction) 
force i S i a. 

Let us proceed to the calculations. In the Lorentz gauge d^A^ = the Maxwell equations take 
the form 

DA^ix) = -Artj^ix) . (14) 

Any solution to these equations may be constructed as the sum of their particular solution and 
a solution to the corresponding homogeneous equations. Although this decomposition is quite 
ambiguous (from the pure mathematical viewpoint) it has a certain physical meaning if one identifies 
the former solution with the own field of the point charged particle, treating the latter as an external 
field describing free electromagnetic waves incident on the particle. The own field of the particle is 
given by the Lienard-Wiechert potentials 



where 



A^'^(x) = -An / G^^\x - z)j,{z)dS , (15) 



G'^'ix) = -^^(x°)(5(x2) (16) 

ZTT 

is the retarded Green's function. Taking together Rels. (H), (|15]), (|l^) we arrive at the standard 
expression for the Lorentz force -F^(s) = F^^{x{s))x''{s) accounting for the self-action of the particle 



F, = + Fi: , Ff{s) = F;:\x{s))x\., , 



s 



s 

F';{s) = = V / 6'{-X'{s,T))x[,{T)X^]{s,T)x'^{s)dT, (17) 

— oo 

X'^(s,r) = xf'is) -x^(r) , 
where F^^^ is the strength of the external electromagnetic field and the square brackets stand for 
the antisymmetrization of indices (without one-half). The abbreviation "rr" labeling the self-action 
part of the Lorentz force points on the common interpretation of this term as that describing the 
radiation reaction, i.e. the back reaction of the radiation emitted by an accelerating charge. 

Note that expression ([17D is not very meaningful as it stands since the integral diverges. The 
divergence arise from the singularity of Green's function (|TBp at the vertex of the future light cone 

= 0, x° > 0. So we need some regularization procedure smoothing the behaviour of Green's 
function G{x) at x = 0. The most simple and efficient way to do this is to replace the 5-function 
entering the expression for G by the following (5-shaped sequence [||: 

S(s) = lim s > 0. (18) 

fi^+o a 

Integrating this expression with a test function on the positive half-line one can obtain the following 
asymptotic expansion: 

-s/a °° 

Sa{s) = '-— = Y,{-ar6^-\s). (19) 



a 

n=0 



Thus, the regularized expression for the radiation reaction force ( [T7| ) can be written as 



(20) 



hm F-(s,a). 



This integral is localized at the point r = as a ^ +0, and hence, can be evaluated by the Laplace 
method . In a small vicinity of the localization point the damping exponential factor behaves as 



,X2(s,s-r)/a 



exp IX X ■ X h 

a a 



1..2 1. - 
-X +-X- X 



+ 0(r^) 



(21) 



Note that for the massive particle < and the asymptotic of the Laplace integral (pO| ) is governed 
by the leading r^-term, while for the massless particle it is determined by r^-term. In the latter 
case 



X 



X ■ X = . 



nr* _1_ /y» , y 'T' 

4 3 12 ' 



(22) 



and the coefficient at is strongly negative if we demand x^ = 1. In such a way we recover the 
positive-curvature condition for the isotropic world lines discussed in the previous section. In the 
massless case the result of integration should have a form of Laurent's series in a^^^ with a finite 
irregular part. The details of calculations are given in Appendix. The result is 



F7(.,a) = -g^ 



Ol .. -3/4 , ^•••2 • -1/2 I 



10 



16 



^ll---2-- I II---. • \ -1/4 

^ 10 ^ 5 ^ 



22 [(5) ,...2... , o--- , /9(^)2 , (5) 2,...2,2\ . 1 , 

-y g i x^+x^x^ + 3x- xx^+ I - X +— X- X +-(x^)^lx^^ + 



(23) 



ar. 



12"/^r (1 + n/4) 



where dots stand for terms vanishing upon removing regular izat ion. The expression (|23|) contains 
three singular terms, two of which are Lagrangian, namely, those at —3/4 and —1/4 powers of a, 
while the rest is not. For example, a~^''^-term is obtained by varying the following functional: 



2 Oil _3 

q —a 4 



X 



-.dr. 



(24) 



As to a~^/'^-term, the corresponding Lagrangian is given by a rather unwieldy expression. In order 
to simplify it let us redefine the Lagrangian multiplier e entering the l.h.s. of the Lorentz equation 
(^ as follows: 

.11 . _ ,. 

(25) 



2 -L-L •••2 -1/4 

e — q a^x a ' 



160 



Then the expression for a ^/'^-term reduces to 



2 "3 (4) _i/4 

q — X u a ' 
^16 ^ 



(26) 



which can be obtained by varying the following functional: 

^2^^-1/4 f (27) 



A practical recipe to restore both the Lagrangians is as follows. We exclude the electromagnetic 
field An{x) from the action functional ([l|) by the Maxwell equations (^). The result is the Fokker 
type |]lOl effective action for the charged massless particle 

Seff = -\l drex' - 2n [ d'x [ dSf{x)G''\x - y)j,{y) = 



{2i 

' dTex^ + '^ ds dr x^'{s)5{X^{s,t))x^{t). 



2 i 2 _ . 

The first term is the usual action of a free massless particle, while the second term describes particle's 
self-action. Using regularization ( |T8|) for the J-function entering the self-action term we may perform 
one of two integrations. The calculations are quite similar to those presented in Appendix. The only 
distinction is that now we cannot use the isotropy condition (^). Instead, we rewrite the exponent 
(pTD in the following way: 



exp 



1..2 1. - 

-x H — X ■ X 



exp x^ X-X— + 0{t^) \ . (29) 

\ a a 



Then the first multiplier is interpreted as a damping exponential factor of the Laplace integral, 
while the second is combined with the rest of the integrand and expanded in the Taylor series 
in r; in so doing, we may keep only terms which are at most linear in x^ and x ■ x, since the 
higher powers will not contribute to the equations of motion 6Seff/Sx = (more precisely, their 
contribution will vanish on the shell of the isotropy condition x"^ = 6Seff/Se = 0). Integrating 
we get Laurent's series in a^^^ without a constant term and with two singular terms giving upon 
variation two of three foregoing singularities. Clearly, the Lagrangians obtained in such a way are 
not uniquely determined as one may add to them any local expression proportional to the squares 
of the Lagrangian constraint i;^ and its differential consequences. 

The reason why performing a self-consistent elimination of the electromagnetic field from the 
quadratic (in A^) action (|T]) we do not reproduce the whole expression for the radiation reaction 
force (p3D including the finite part and the a~^/^-term is as follows: Only symmetric part 

G-*(a; -y) + G^^\y - x) = -5{{x - yf) (30) 

of the retarded Green's function (p!6| ) actually contributes to the effective action (^), while sub- 
stituting the Lienard-Wiechert potentials to the Lorentz equation (|) we use, in fact, the entire 
Green's function. The symmetric and antisymmetric parts of the Green function lead to the dif- 
ferent groups of terms in (^); the first terms are invariant under reversion of particle's trajectory 
T — r, whereas the second ones acquire the minus sign. Therefore only the invariant terms can 
arise from the effective action (|28|) . 



According to general prescriptions of the renormalization theory the singular coefficients a 

the regularized expression (^) are replaced by finite constants to be fixed from an 
experiment. Since there is no evidence for a preferable choice for these constants we may put them 
to zero to simplify further consideration. Then the radiation reaction force FJ^"^ is given by the 
second line in Eq.(^). The non-Lagrangian nature of this force manifests itself in the violation 
of the time-reversion symmetry, r — r, which is consistent with irreversibility of the radiation 
process. 

In conclusion of this section let us note that although the expression (|23| ) for the radiation reac- 
tion force was derived in a special gauge, it is not hard to rewrite it in an arbitrary parametrization. 
For this end one should just treat the overdot derivatives as invariant ones, 

d = ^:j-, (31) 



where r is already an arbitrary evolution parameter. Note that in an arbitrary parametrization 
the effective equations of motion for a charged massless particle involve a sixth derivative (for 
comparison, the Lorentz-Dirac equation for d=4 massive particle is of the third order). 



4 Reduction of order ^ 



Let us write down the effective equation of motion governing the dynamics of a massless charged 
particle, 

D{iDx>') = n^Dx"" - — (D^x^ + {D^xfD^x^" + 3D^x ■ D^xD^x^'+ 

5 



+ 



^{D'xf + —D'x ■ D'x + -{iD''x ff 
7 7 5 



(32) 

Dx^' ] , {Dxf = , 



where e(r) = q~^e{T)\^ is a nonzero scalar function and = ri^^F^^*{x). 

Since the order of the equation is greater than two, it cannot be assigned with a straightforward 
mechanical interpretation: In the realm of Newtonian mechanics a state of the particle is unambigu- 
ously determined by specifying its position and velocity, that is obviously insufficient to extract a 
particular solution to our equation. This problem is analogous to that of a mechanical interpretation 
of the Lorentz-Dirac equation which also has too many solutions, but not all of them are physi- 
cally meaningful. To extract the subspace of physical solutions we impose an additional selection 
rule 0]: The physical solutions to the Eq. (|3^) are only those which have a smooth limit upon 



switching off the interaction, i.e. when q 0. In other words, we treat the radiation reaction as a 

small perturbation deforming the conventional dynamics of the charged particle (i.e. the dynamics 

without account of the self-action) rather than something introducing extra degrees of freedom to 

the theory. This treatment is justified by a small parameter q at the higher-derivative terms in the 

^We are thankful to S.L. Lyakhovich for illuminating discussions on various aspects of the reduction-of-order 
procedure for higher-derivative Lagrangian systems. 



r.h.s. of the Eg. (P^ . Instead of seeking for solutions which can be smoothly continued to g = 0, 
one may equivalently seek for a second-order differential equation (with a smooth dependence of 
q) any solution of which would be a solution to the initial equation. The general procedure for 
obtaining such an equation is known as the reduction of order. It consists in successive elimination 
of the higher derivatives from the equation (|3^) with the help of all its differential consequences; in 
so doing, the order of the equation increases at each step of the procedure, but the order of q at 
higher derivatives increases as well, so that, in the limit, we get a power series in q with coefficients 
depending on x, Dx and D'^x. The implementation of this procedure to the Lorentz-Dirac equation 
may be found, for example, in [|], |^, ||. It should be noted, however, that the resulting equation 
for the massless particle appears to be unresolved with respect to the highest (second) derivative 
that poses the questions about existence and uniqueness of its solution given the initial position 
and velocity of the particle. Without going into detail we just note that for a homogeneous external 
field {Q = const) the reduced equation can be perturbatively resolved with respect to D'^x with the 
right hand side given by a power series in q. Up to the first order in q it reads PI 



x"" 



2q {n^xY xQ^x.^^..^ 
5 V e^ ^ ' 



-e-^{^lxY+ 



_ Qxn^x lixn'^xf 

+ h --^ — 



(33) 



0{q' 



where overdot stands for the invariant derivative D and the function e, being determined from the 
identity i;^ = 1, is given by 

-xn'^x + 0(g2) . (34) 



As above we are restricted to initial velocities obeying condition e 7^ (compare with ([T0|) ) which 
ensures existence of a solution for a sufficiently small time interval. In the gauge e = 1 the above 
equation takes the form 



x^" = q{nxY+ 



+- 



2q^ f{n^x) 



xVL'^x 



xVL'^x 
[iVfix)'' 



{n^x) 



9 xVL^x lixVL^x) 



(35) 



+ o{q') 



where overdot denotes now the ordinary derivative w.r.t. the evolution parameter. Note that 
although the condition {xVl^x) < has been assumed at each stage of our derivation the resulting 
equation has the correct uniform limit upon switching off the external electromagnetic field (i.e. if 
we set Q = eQ' and then let e ^ 0). The last statement remains true even with account of higher 
orders in q, so the limiting equation x'^ = describes the free motion. 

Thus, the perturbative treatment of the radiation reaction, being apphed to the case at hands, 
leads to the meaningful second order equation at least for the constant fields. 



^Here we use matrix notation: (ili)^ = fl'^x'^ , xfl^x — ii.fl^rt'^x'^ , and so on. 



n 



5 Conclusion 

To summaries, in this paper we have derived the effective equations of motion for a massless charged 
particle interacting with external electromagnetic field as well as its own one. Although the approach 
we follow here is quite similar to that used in the case of a massive particle 0] the final results 
considerably differ from each other. First, not all the divergencies are Lagrangian, i.e. may be 
canceled out by adding appropriate counterterms to the original action functional; even though such 
counterterms do exist (for some divergences) they are not uniquely defined. The ambiguity reflects 
the freedom to add to the Lagrangian any local expression proportional to the square of the isotropy 
condition and its differential consequences. The non-Lagrangian divergences in turn are recognized 
as those breaking up the time-reversion symmetry. Second, consistency of the dynamics imposes 
some restrictions on initial data of the particle and/or the form of external fields. Third, after 
renormalization we get a fifth order differential equation, instead of third-order equation derived by 
Dirac for a massive particle. 

All these distinctions have a common origin related to the special local structure of the isotropic 
world lines: The interval between two nearby points x{t) and x{t + 6t) on an isotropic curve, lying 
in a general position, is proportional to 5r^, rather than to be linear in 6t as is case for the time-like 
world lines of massive particles. This leads to the stronger singularities, as compared to the massive 
case, and as result to a different structure of the radiation reaction force. 

Finally, the straightforward application of the reduction-of-order procedure leads to the differ- 
ential equation in a form unresolved with respect to the second derivative which may lead to some 
difficulties upon its integration. Besides, it may have no well-defined limit upon switching off an 
external electromagnetic field as seen from the structure of the first radiative correction (^). The 
last fact may indicate a certain instability of the classical dynamics and deserve a further study. 

It would be interesting to re-derive our equation on the basis of usual energy conservation 
arguments in order to gain a more physical insight into the problem as well as an independent test 
for consistency of the formal renormalization procedure we have used. 

Appendix 

The regularized expression for the self-action force (^OD involves the following integral: 



Here we apply the Laplace method to obtain the Laurent series for the function with respect to 
a^/^. More precisely, we are interested only in constant and irregular parts of this series (contribut- 
ing, respectively, to the finite and the divergent parts the self-action force) since the regular terms 
vanish upon switching off the regularization. 

In order to simplify formulae we use standard notation from the linear algebra for the inner and 
the exterior (skew-symmetric) products of two vectors: a ■ b = a^^h^, a Ab = {a^h^ — a^h^). Thus, 



oo 




(36) 







1 n 



the pre-exponential factor in the integral ( P^D may be written as 

oo 

x{s - r) A X{s, s-t) = J2 r'^bUs) , (37) 



where 



In particular, 



m=2 



m (") (m-n+1) 
X A X 



^ n\(m — n)\ 

n=l ^ ' 



1 (2) (1) , 1 (3) (1) 

62 = --a;Aa;, 63 = -xAx, 

1 (4) (1) 1 (3) (2) , 1 (5) (1) 1 (4) (2) 

Oa = — X A X X A X , = — X A X H x A x , 

8 12 '30 24 



1 (6) (1) 1 (5) (2) 1 (4) (3) , 1 (7) (1) 1 (6) (2) 1 (5) (3) 

Ofi = xAa; x A x x A x , 67 = x A x H x A x H x A x . 

^ 144 80 144 ' ^ 840 360 360 

Using the basic identities 

x2 = , X^ = 1 

as well as their differential consequences 

(4) 

x-x = 0, x-x = 0, x-x = — 1, x-x = 0, 

(4) (2) (3) (5) (1) (3) (5) (2) (4) (3) (6) (1) (4) (3) 

/■>'■ /■>'■ _[_ /-y» ^ /-y» /■>'■ /-v» ^ /-v» /y» _[_ J /-v» rr* 'If 'If ^ If If I I 

the argument of the exponent (|36|) can be written as 

oo 

s - r) = (x(s - r) - x{s)f = ^ c^r^ = 



where 



fe=4 

1 (2) (2) 4 1 (3) (3) g 

= X ■ X r H X ■ X r — 38 

12 360 ^ ^ 

1 (3) (4) / 1 (3) (5) 1 (4) {4)\ 8,8. 

X ■ X r + X ■ X H X ■ X r + o r , 

360 V 1260 1344 ' ^ ^' 



fc-l (n) (fc-n) 



^ — n)! 

Substituting ( |37D and (|38| ) to the integral (|36| ) and making replacement t ^ a^/^t we find that 



^^^""^ ^ ^'^Ya j exp a("-^)/V"c„ | + (regular terms) (39) 



The expansion for the exponent reads 

e-Vi2(i + ^1/2^6^^ ^ ^3/4^7^^ ^ ^[^8^^ ^ T^^cll2]) + o{a) . (40) 

Muhiplying this expression on the pre-exponential factor 



7 

m — 3 



m=3 

and contracting with i;'' we get 

+a[T'^h + (cefes + C764 + cgfea)!-^^ + Ic^fegri^])^^ + o(a) . 
Substituting the last expression into (|36|) and using the integration formula 



00 



e~^Vi2^n^^ ^ 12^r { ) (41) 



4 



we finally arrive at ([23|). It is interesting to note that the resulting expression for F"'{s,a) turns 
out to be orthogonal to both and x'^. 
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